Spin-orbit induced mixed-parity pairing in S^RuO^ a quantum many-body 

calculation 
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The unusual superconducting state in S^RuCU has long been viewed as being analogous to a 
superfluid state in liquid 3 He. Nevertheless, calculations based on this odd-parity state are presently 
unable to completely reconcile the properties of Sr2RuC>4. Using a self-consistent quantum many- 
body scheme that employs realistic parameters, we are able to model several signature properties of 
the normal and superconducting states of S^RuCU. We find that the dominant component of the 
model superconducting state is of even parity and closely related to superconducting state for the 
high-Tc cuprates although a smaller odd-parity component is induced by spin-orbit coupling. This 
mixed pairing state gives a more complete representation of the complex phenomena measured in 
SnRuCU. 



The unusual electronic properties of the layered ma- 
terial Si - 2Ru04 remain of significant interest eight years 
after Mackenzie and Maeno's review [1] of an already 
extensive literature. Experiments show a quasi two- 
dimensional Fermi liquid state for T < 100K with an 
unconventional superconducting state emerging at T c = 
1.5K. A qualitative description of the normal state Fermi 
liquid has emerged through band structure calculations. 
Three quasiparticle bands, a, /3 and 7, originating from 
d X z, d yz and d xy orbitals of tetragonally-coordinated Ru 
ions, cross the Fermi surface [2, 3]. The nature of the su- 
perconducting state is less well established, but the lack 
of in-plane Knight shifts for T < T c [4] and the pres- 
ence of finite signals in symmetry sensitive probes such 
as the polar Kerr [5] or the Josephson effects [6, 7] point 
to spin-triplet, odd-parity pairing and, more specifically, 
to a chiral pairing state with zero total pair spin per- 
pendicular to the planes, i.e S z — 0. Nonetheless, other 
experiments are not so well-explained by such a pairing 
state [8]: the specific heat for T << T c [9] shows evidence 
of nodes or minima in the superconducting gap function, 
but a chiral p-wave state is nodeless; the out-of-plane 
Knight shift [10] is weakly temperature dependent be- 
low T c ; and edge currents are not detectable at a level 
that is two orders of magnitude below the expected value 
[11]. More recently, signatures of half-quantum vortices 
below T c were explained in terms of equal spin pairing, 
i.e. S z / [12]. Thus, the nature of the bulk super- 
conducting state is still not resolved completely and it 
is essential to test current working assumptions with mi- 
croscopic calculations that are as detailed as possible. 

Our calculation starts with a tight-binding representa- 
tion of an LDA band-structure for Sr 2 Ru04 as developed 
by Pavarini and Mazin [13]. Low energy electronic states 
are represented with a basis consisting of a d xy , d xz and 
d yz orbital for each Ru ion and these are 2/3 filled on av- 
erage. Ru ions are strongly coupled within planes where 
they form a square lattice of side a = 3.87A. Weak in- 
terplanar couplings corresponding to the z direction are 
ignored. For the Ru atomic spin-orbit interaction, A s ■ I, 



we use A — 70 meV [14] . The interaction of atomic spin 
and orbital moments and a uniform external field is given 
by -ptB(g s s + l) ■ B. 

Dynamical correlation effects are modeled on the ba- 
sis of an atomically local bare interaction with three pa- 
rameters; U and U' for inter- and intra-orbital Coulomb 
repulsion and J for intcrorbital exchange; and the fluctu- 
ation exchange approximation (FLEX) . Here we will only 
consider fixed ratios of U' /U = 2/3 and J/U = 1/6 that 
are typically used in these calculations [15]; variations in 
these ratios have little impact on general features of our 
results. Having fixed the ratios we find that U < 1 eV is 
reasonably consistent with the observed T c — 1.5K, but 
we will in some instances use a somewhat larger value of 
U = 1.2 eV as the higher T c that results makes super- 
conducting state calculations more feasible. To make the 
scheme computationally feasible, the dynamical cluster 
approximation [16] is used to reduce dynamical correla- 
tions to a range of L c << L where L ~ 100a is the physi- 
cal lattice dimension and we numercaly demonstrate that 
our results well-describe the large L c limit. Translation- 
ally invariant, but otherwise random initial configura- 
tions are used to seed the self-consistent procedure which 
then terminates when a stable state is produced. 

Figure 1 compares the renormalized normal state (T = 
60K) electronic structure of S^RuCU as determined by 
FLEX calculations for U = 0.9eV with bare band disper- 
sions derived from LDA [13]. The FLEX band disper- 
sions have significantly smaller slopes at the Fermi level 
than unrenormalized LDA counterparts. These indicate 
a significant mass enhancement for each band that is ex- 
pected from specific heat [17] and dHvA [18] measure- 
ments. The imaginary part of the self energy for each 
band is proportional to (E — ef) 2 as is expected for a 
Fermi liquid. 

As the calculations can be used to determine the elec- 
tronic structure with both energy and momentum res- 
olution, it is natural to present a comparison to angle 
resolved photoemission experiments. Fortunately, once 
the problem of a surface reconstruction was resolved [19], 
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FIG. 1. (Color online) Comparison of low energy band struc- 
ture derived from FLEX calculations (symbols) and LDA 
(solid lines) along the indicated symmetry directions. The 
vertical lines through each point correspond to the inverse 
lifetime of the state. There is a significant renormalization of 
each band, with the 7 band showing the strongest renormal- 
ization along T-M. 



function, a \T — a2| ai to estimate T c = 02. Results for 
T c versus U, A and L c are shown in Table 1. While 
the numerical parameter L c does affect quantitative re- 
sults, qualitative features are reasonably well described 
with L c — 4a. With spin-orbit coupling (A = 0.07) we 
find that T c ~ 20 - 30K for 17=1.2 eV and T c < 10K 
for U = 0.9 eV. Thus, U < 0.9 eV is most consistent 
with the observed T c , unless the spin-orbit interaction 
for Sr2Ru04 in actuality is much larger than 0.07 eV. 
Recall that [7=0.9 eV leads to quantitative agreement 
with photoemission results as well. 
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there have been a plethora of high quality photoemis- 
sion experiments that probed the bulk normal state elec- 
tronic structure of Sr 2 Ru04 [20-23] . For the f3 and 7 band 
crossings along F-M, the calculated dynamical renormal- 
izations of the Fermi velocity, vf,lda/vf,flex where 
vp = dE/dk\k=k F , are 1.9 and 3.4 which compare well to 
the experimental values of (approximately) 1.7 and 3.3 
[22]. 

A more striking renormalization is the shift of the 
van Hove singularity in the 7-band at the M point from 
90 meV to 20 meV in going from the LDA to FLEX re- 
sults. This feature is clearly apparent in photoemission 
experiments utilizing either thermally occupied states or 
doping [24]. This confirms that FLEX is capturing non- 
trivial normal-state correlation effects that could be key 
precursors to the superconducting transition. We note 
that spin-orbit coupling leads to non-negligible spin and 
orbital mixing for 7-band quasiparticles at tp along T- 
M with that d xz orbitals contributing about 15% of the 
quasiparticle weight. Since the spin-orbit induced mix- 
ing of d xy and d xz orbitals involves a spin flip, i.e. 7 + ,s~, 
these quasiparticles have a wavefunctions of the form 
\1p1p >^ ci\d xy , J \> +c-2\d xz ,\.>, i.e. they are not pure 
spin-eigenstates. 

Superconductivity is associated with the stabil- 
ity of finite anomalous Green's function amplitudes, 
F v<tuI(t ,{R,t) = -% < <& a {$.,T) 4, CT ,(0*,0) > . As an or- 
der parameter amplitude for the superconducting state 
we take m 2 p (T) = Yj VaV 'a'R \F V aWa'{R,T = 0)| 2 and fit 
the temperature dependence of m p (T) to a power-law 



TABLE I. The transition temperature, T c , and pair wave 
function weights by spin, p( s ' s ^) i as a function of the 
Coulomb repulsion strength ([/), and spin-orbit coupling con- 
stant (A) and dynamical cluster size (L c ). In the absence of 
spin-orbit coupling (A = 0), the pairing state is a pure singlet. 
Spin-orbit coupling induces triplet components to the pairing 
state, primarily between electrons with equal z-components 
of spin, S z — ±1. 

To examine the spatial, orbital and spin dependence 
of the pairing state, we define a normalized pair wave 
function and decompose it into eigenfunctions of total 
spin, (S, S z ), 
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from which we obtain a fractional weight by total spin, 
p(S,S z ) ^ ky taking the complex squares of ^lj, (R) and 
tracing over orbital indices and displacements. In Table 1 
we show results for these weights as a function of U, A 
and L c evaluated at T ~ T c . These results show a pure 
singlet pairing state when spin-orbit coupling is absent, 
i.e. A = 0. Further examination of the orbital and spatial 
properties of ^! v „i (R) show that this singlet state has 
d x 2_ y 2 symmetry with respect to rotations in the plane, 

i.e. ^ X y-J y {ay) = —Skiy-J y (ax) , and that the only sig- 
nificant orbital contributions come from the d X j,-orbitals 
with the very small contributions from d xz and d yz or- 
bitals vanishing completely in the limit J — > 0. Conse- 
quently, for A = we find that the superconducting state 



leaves the a and /3 bands ungapped or very nearly so 
and the 7-band gap has nodes typical of a d x 2_ y 2 pairing 
state. 

A finite spin-orbit coupling constant, A ^ 0, finite 
triplet pairing components. The largest amplitudes for 
triplet pairing occur for equal spin pairs, i.e. S z = ±1. 
Further examination of the wavef unctions for S z = ±1 
reveals that they describe interorbital pairing between 
an electron in the d xy orbital and an electron in cither 
of the d xz of d yz orbitals. Also we find ±7r/2 phase 
shifts between the two types of interorbital terms in this 
wavefunction with opposite relative signs for the up and 
down spin pairs, i.e &xy,yz(ay) = s~ m ' 2 ^ X y,xz(ax) and 
^xy^yl (ay) — e +llT ^ 2 ^ xy~Jz (ax). Odd-parity states of 
this form are "topological," time-reversal invariant su- 
perconducting states and their unique properties have 
recently been described by Qi, et al. [25]. Tada et al. 
[26] have performed model calculations that point to the 
possibility of generating such a state at the interface via 
spin-orbit coupling although they suggest that weaker 
spin-orbit effects lead to a chiral p-wave state in the bulk. 
Note that for such a state up and down spin equal-spin 
pairs have opposite chirality, so no net charge current is 
expected on the edges, a result that is consistent with 
experiment [11]. We note that for fixed values of U and 
A, the relative contribution of singlet and triplet pairing 
terms is weakly dependent on the numerical parameter 
L c with little change observed in going from L c /a = 8 to 
16. 

In Figure 2, we show the quasiparticle excitation spec- 
trum for T ~ 0.5T C along the T-M cut evaluated with 
U = 1.2 eV, A = 0.07 eV and L c = 8a. The charac- 
teristic excitations on both sides of the superconducting 
gap are evident in the 7-band near k x — n/a. No ex- 
tra gap is readily apparent in either the a or /3 bands, 
or in any of the bands along the zone diagonal. This 
is the structure observed for pure d-wave pairing when 
A = and while new, but smaller gap features are ex- 
pected with the triplet component, they are not readily 
apparent given the momentum/energy resolution (as in- 
dicated by the plotted points) in this calculation. Recall 
that nodal structure in the superconducting gap observed 
in these FLEX results is consistent with observations for 
the specific heat superconducting Sr 2 Ru04 [9] whereas a 
chiral p-wave state is expected to be nodeless. However, 
it remains to be demonstrated that a superconducting 
gap function of the form shown here is consistent with 
all probes that are sensitive to the quasiparticle excita- 
tion spectrum below T c , all of which are consistent with 
a nodal strucutre, but as of yet, of undetermined form. 

Next we consider the temperature dependence of the 
Pauli spin susceptibility. Measurements of the Knight 
shift for < T < T c suggest a weak variation with tem- 
perature for both Xxx and Xzz-, the in-plane and out- 
of-plane components of the uniform Pauli susceptibility 
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FIG. 2. Quasiparticle spectrum along the T-M direction for 
T < T c . The usual signature of a superconducting gap ap- 
pears in the 7-band. There is no evidence for a smaller gap 
forming in the /3-band along this cut to within the momentum 
resolution of the calculation. 



[4, 10]. Calculated susceptibility results for U = 1.2 eV 
and L c = 4a as a function of reduced temperature, 
T/T c , are displayed in Figure 3. For A = we find 
that Xxx — Xzz is suppressed by about 1/3 below T c , 
but the curves do not approach zero as T — > as the 
a-band and /3-band Fermi surfaces remain ungapped. 
In the mixed pairing state, for A = 0.07 eV, Xxx and 
Xzz are no longer equal, but each has a weak tempera- 
ture dependence below T c with changes of less than 10% 
compared to T ~ T c . Thus, the qualitative behavior 
of the temperature-dependent spin susceptibility for this 
mixed pairing state apparently is in good agreement with 
Knight shift data. Further, these results suggest that the 
small Knight shifts in Si'2Ru04 are caused by the mixed- 
spin nature of quasiparticle states, a general mechanism 
for small Knight shifts described by Anderson [27], and 
not by a bulk p-wave symmetry pairing state. 
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FIG. 3. Temperature dependence of the Pauli spin suscepti- 
bility, Xsp, for zero and finite spin-orbit coupling constant, A. 
For A = the equal in-plane and out-of-plane components are 
suppressed by the superconducting transition In contrast, for 
A = 70 meV, both the in-plane and out-of-plane components 
are suppressed by less than 10% for T < T c . 

We note, though, that the overall scale for the calcu- 



lated susceptibilities are a factor of about four too small 
in comparison to experiment [17]. The lack of quantita- 
tive accuracy in this instance is most likely a failure of 
FLEX as the approximation has been shown to underes- 
timate the large enhancements of the magnetic suscep- 
tibility that are expected in the vicinity of a magnetic 
phase transition for the special case of the single-band 
Hubbard model [28]. 

Thus, the mixed pairing superconducting state is con- 
sistent with major experimental features of Sr 2 Ru04 
in both the superconducting and normal states. How- 
ever, the time-reversal symmetry preserving nature of 
this mixed-parity state places these results at odds with 
the interpretation of polar Kerr and tunneling experi- 
ments [5, 7]. We note that such an interpretation of the 
tunneling data requires the presence of superconducting 
domains of size < 1/im, which is too small to account for 
a finite signal in the infrared frequencies used in the polar 
Kerr effect experiments [8] . Also the existence of a finite 
polar Kerr effect at the pseudogap temperature, T*, in 
the cuprate superconductors suggests that this probe is 
sensitive to correlations other than those directly related 
to a time-reversal symmetry breaking superconducting 
state [29]. It remains to be seen whether the correlation 
between experiment and calculations can be improved by 
accounting for additional details such as vortices, defects, 
impurities, surface effects, domains and other perturba- 
tions, all of which have been found to be necessary for 
enabling the theory for a pure p-wave superconducting 
ground state to be consistent with the wide range of ex- 
perimental studies of the system. 

The results shown here are supported by earlier the- 
oretical work. For example, Kuwabara and Ogata [30] 
showed that fluctuations may favor a d x 2_ y 2 pairing given 
certain parameters in an effective one-band model for 
Sr2Ru04. In the course of developing a theory for a 
p-wave superconducting state for Sr2Ru04, it has long 
been understood that spin-orbit coupling might play a 
role in stabilizing one of the possible p-wave pairing states 
[14], but less attention has been placed on the impact of 
spin-orbit coupling on a nominally d x 2_ y 2 pairing state. 
However, at the conclusion of their study of spin-orbit 
effects on the band structure of Sr 2 Ru0 4 and Sr 2 Rh0 4 , 
Haver kort, et al. state "In fact, singlet and triplet states 
could be mixed, blurring the distinction between spin- 
singlet and spin triplet-pairing ... [31]." Indeed, we find 
that when dynamical correlations and spin-orbit inter- 
actions are simultaneously taking into account using re- 
alistic parameters for Sr 2 Ru04 such a mixed-parity bulk 
state is realized. As FLEX is an approximate many-body 
scheme its identification of the superconducting state for 
the underlying three-band model must be considered ten- 
tative until it is computationally feasible to include terms 
that are missing in FLEX (i.e. vertex corrections in the 
self-energy series) or perform a quantum Monte Carlo 
analysis. Nonetheless, this FLEX analysis represents is a 



methodological step forward as it incorporates physical 
effects not included in earlier finite-order perturbative 
studies [32, 33], namley (i) collective spin fluctuations 
via an infinite set of ladder diagrams and (ii) dynamical 
quasiparticle renormalizations via self-consistency of the 
electron Greens function. As these effects are not small, 
it is unsurprising that a different conclusion results with 
respect to the pairing symmetry. 
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